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Abstract 

We present details of our evaluation of the NNLO QCD corrections to B{B — )■ Xs'^) in the 
heavy charm limit [rrtc ^ ^). Results of this calculation have been essential for estimating 
O (a^) effects in this branching ratio via interpolation in mc- 



1 Introduction 

The inclusive branching ratio B{B — )■ Xg'j) is well-known to provide important constraints on 
extensions of the SM [1] . Its evaluation is based on an approximate equality of the hadronic 
and and partonic decay widthq^ 

r(5 ^ x,7),^>,^ ^ r(6^xf7),^,,^, (1.1) 

where X^ stands for s, sg, sgg, sqq, etc. This approximation works well only in a certain range 
of Eq, namely when Eq is large {Eq ~ ^) but not too close to the endpoint {mb — 2Eo S> Aqcd)- 
It has become customary to use Eq = 1.6 GeV ~ ^ for comparing theory with experiment. 
Calculations including the C(af) and C(aem) effects in the SM give |11[5] 

BiB ^ X,7),,>,ec.v = (3-15 ± 0.23) x 10-^ (1.2) 

where the uncertainty is dominated by O^asAQco/nif,) non-perturbative effects [2]. 
The currently available experimental world averages read 

r.(r.^y . ^ / (3-55 ± 0.24exp ± 0.09^odci) X 10-4 |6], 

^ ^«^^-.>i.aGcv I ^3_5Q ^ Q_^4^^^ ^ O.lO^odci) X 10-4 m- ^ ' 

They have been obtained from the measurements of CLEO [8], BABAR [9] and BELLE [TO|[TT| 
by extrapolation in i^n according to various photon energy spectrum models, whose parameters 
have been fit to dataci The SM prediction (II. 2p and the averages ( II. 3p are consistent at the 
1.2a level. 

The O(a^) contributions to the branching ratio amount to around 10%, which exceeds the 
experimental errors and theoretical non-perturbative uncertainties. However, these corrections 
have not been included in a complete manner in Eq. (II. 2p because their charm-mass dependence 
remains unknown beyond the BLM-approximation [12]. Instead, we have calculated all the Trie- 
dependent non-BLM corrections in the rric ^ ^ limit, and performed their interpolation in 
rUc down to the measured value rric — ^, assuming that they vanish at rric = 0. Our previous 
paper [5] contains only the final analytic expressions for the large-mc results together with a 
description of the interpolation. Presenting details of the large-mc calculation is the purpose 
of the present article. 

The paper is organized as follows. Sec. [2] is devoted to recalling the relevant definitions 
from Ref [S]. Sec. [3] contains an explanation why we did not use asymptotic expansions of 
three-loop on-shell Feynman diagrams. Our actual method that involved charm decoupling at 
the Lagrangian level is described in Sees. H] and [51 We conclude in Sec. [61 Expressions for the 
relevant functions (plj that originate from b — )■ s'jg are collected in the Appendix. 



^ Corrections to Eq. (|l.ip of various origin have been widely discussed in the literature, most recently in 
Ref. [2] ■ A compact overview of previous results can be found in Ref. ^ . 

^ Ref. [6] gives a larger error than [7] because it includes results at Eq > 1.8 GeV from the older measure- 
ments [Slini only, ignoring the more precise ones from Ref. pTj . 



2 Notation 

We shall strictly follow the notation of Ref. [5]. The present section collects the most important 
definitions only. The effective Lagrangian that matters for evaluating QCD corrections to 
b^XP-f readfi 

AG * 

^cs = >Cqcdxqed(^,c?, s,c,6) H y^V*,VtbJ2CiQi^ (2-1) 

where the local flavour-changing operators Qi arise from decoupling of the W boson and all the 
heavier particles. We shall need explicit expressions for the following ones: 

Ql = {sLltJ''CL){cLYT''hL), Q2 = {sLlf,CL){cLrbL), 

Ql = jt^nn,{sLa^^-hR)F^,, Qs = ji^mk{sLa^''T%R)G%. (2.2) 

The remaining three (Q3, Q^ and Qe) are similar to (^4 but involve different Dirac and color 
structures. The sum over q in (^4 runs over all the active flavours (u, d, s, c, h). Masses of the 
light quarks (u, d, s) are neglected throughout the paper, except for the coUinear logarithm in 
Eq. (lA.ap of the Appendix. 

The Wilson Coefficients (WCs) Ci are assumed to be MS-renormalized at the scale //(, ~ ^. 
To avoid scheme- dependence at the Leading Order (LO) in QCD, one usually works with certain 
linear combinations called "effective coefficients" 

{Ci, for z = 1,...,6, 

Ci + Y^Uvfi,, forz = 7, (2.3) 

where yj = (0, 0, — |, — |, — y , — ^)j and Zj = (0, 0, 1, — |, 20, — y)j in dimensional regular- 
ization with fully anticommuting 75. In the SM, all the Cf^{fib) are known up to 0{a'j.) [T3HT5] . 
We are interested in evaluating the Next-to-Next-to-Leading Order (NNLO) QCD corrections 
to the ratio of partonic radiative and charmless semileptonic decay rates 



T[b ^ x,^]e,>eo v:ytb ' 
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6«.n. ^ ^e.^^^) cf{^^,) K,„ (2.4) 



^ ^,j=l 



where the symmetric matrix Kij is perturbatively expanded as 

Kij = S„5,r + a.K^f + a^.Kg^ + O [al) with a, = ^^^^. (2.5) 

Splitting K^j into the BLM {K^j °) and non-BLM (K^j '^) pieces is done in a standard 



manner: 



'(2) _ /I ^ I R _ 7>'{2)/3o I T^{2)vem 



Ki;> = A.uf + B,, = K^;>'^ + K^;>'-, (2.6) 

^ Including the electroweak or CKM-supprcssed corrections would require adding more terms to Eq. (|2.ip . 



Figure 1: b ^ s'yqq {q = u, d, s) diagrams with (5i,2 vertices that survive in the large-rric limit. 



where rif stands for the number of quark flavours in the effective theory (12. ip . and 
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(2.7) 



In the 0(5^) correction, one can safely ignore the small C3 4 5g(/ib). Thus, it is sufficient to 

with such indices are concerned, all of 

A calculation 



consider K^j with i,j E {1, 2, 7, 8} only. As far as Kl 



them except i^jg 



K2)/3o 



(2)/3o 
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and -^^28 ^^6 known [T6HT9] for the measured value of m. 



of -ft'is and K2Q is underway 

Effects related to the absence of real cc production in 6 — )■ Xf 7 and to non-zero masses of 
b and c quarks in loops on gluon propagators belong to i^^-j . They are presently known 
for all the i,j G {1,2,7,8}, either for arbitrary rric or at least in the vicinity of its measured 
value [ITlEl]. In fact, charm quark loops on gluon propagators are the only source of nic 



'(2) 7^(2) 



'(2) 



dependent terms in -ft'77 , K^s and Kgg . Therefore, the only quantities for which the Trie- 
interpolation still needs to be performed are i^jj ^^™ for i G {1,2} and j G {1,2, 7, 8}. In the 
following, we shall restrict our considerations to those cases only. 

Before closing this section, let us remark that our large-?7ic calculation is not 100% complete. 
There exist certain simple though yet uncalculated contributions to ii'|- °™ with i G {1, 2} and 
j G {1,2,7,8} that survive in the large-?7ic limit. They originate from the four diagrams in 
Fig. [T]that may interfere either with b — ?■ s'-fqq contributions of Qt.s or just with themselves. 
Their effect on the decay rate is of order 5^, and it is expected to be numerically very small 
due to limited four-body phase space left out by the high photon energy cutoff -Eq ~ ^- A 
convention advocated in Ref. [TS], which we follow here, is to exclude those uncalculated terms 
from the BLM contribution, even though some of them are proportional to the number of 
massless flavours. We shall comment on this issue again in Sec. |5l 



3 Choice of the method 



^(2)rcm 



for i G {1, 2} and j G {1, 2, 7, 8} in the m^ > ^ limit. 



Our goal amounts to evaluating K^j 

at the leading order in — |. On general grounds, one expects results of the form 
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Figure 2: One of the on-shell non-planar diagrams contributing to ^27 



where X^- are ttIc- independent. A straightforward method to perform such a computation 
via asymptotic expansions [22] would involve three-loop on-shell vertex diagrams like the one 
shown in Fig. [2] in the heavy-charm limit. Application of the so-called hard-mass procedure 
to such diagrams leads to one-, two- and three- loop vacuum integrals with mass scale rric, as 
well as one- and two- loop on-shell vertex integrals with external momenta pi = ml, pf = 
and p'i = 0. Internal lines of the vertex diagrams are either massless or carry mass m^. At the 
two-loop order many different cases occur since up to four bottom quark lines can be present, 
and the photon can couple in all possible ways to the charm or bottom quark. In 2006, i.e. at 
the time when our actual calculation |5] was performed, some of the relevant two-loop on-shell 
massive vertex integrals remained unknown. Furthermore, in addition to the virtual corrections, 
also contributions from real gluon radiation had to be considered, which involved phase-space 
integrals in parallel to the loop ones. 

Although technical challenges related to the asymptotic expansion method are certainly 
manageable, we decided to follow a field theory based approach, as already mentioned in the 
Introduction. This method takes advantage of the fact that charm decoupling at the Lagrangian 
level can be facilitated with the help of Equations of Motion (EOM). In effect, all the necessary 
two-loop on-shell vertex integrals could be reduced to the (planar) ones that are already known 
from the Next-to-Leading Order (NLO) calculations of 6 — )■ S7 [23] . 

Before discussing in more detail the charm quark decoupling in the next section, let us recall 
the large-TTic results for K^j with i G {1, 2} and j G {1, 2, 7, 8}. Contributions to them from 
the two-body channel b ^ s'-f have been evaluated with the help of asymptotic expansions using 
the program exp [2UI2S]. Vacuum integrals were treated with MATAD [2B], and the reduction of 
two-loop vertex contributions to master integrals, which can, e.g., be found in Ref. [27], was 
performed with the help of AIR [28]. AIR is written in MAPLE and is based on the Laporta 
algorithm [29]. A more flexible and more powerful alternative, which is available since 2008, 
would be the program FIRE [50] . 

Following Ref. [3], we write 

Kir" = -QKir^ = (3o Re {-lr^^\z) + 2 [a{z) + b{z) - f ] L, - f^L^} + 24>T\S), (3.2) 
Kif^' = 2 (1 + 5.,) 0Sf °(5), for z, J ^ 7, (3.3) 



where Lh 
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expansion of Re rg 



(2), 



and 6=1 — =^. Our results for the first three terms of the large- 2; 
z), Rea{z) and Re 6(2) read 



Rerf^' 



8 1^2 2 + 112 In 2 + 27650 
g lii ^ -r 243 iii ^ T^ 6561 



if^ln^^ 



2 V405 



572 
18225 



In^; + 



10427 
30375 



135" 



+ J_ f_86_ 1^2 z 



22 ^2835 



1628 
893025 



Inz + 



19899293 



125023500 



-405^ 



om 



Reaf^;) 



Re h{z) 



3 in xJ -^ g -^ ^ 1 27 "^ ^ ^ 486 J 



+ ^(l^ln 



2; + 



1393 
24300 



O 



_Ainr + A _ 1 



81 



81 



45 ^^^ ^ ^ 2025 y 
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487 
33075 
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(3.4) 



which has been confirmed in Ref . [17] using a numerical evaluation of the coefficients at z ^ In" z 
{k,n = 0, 1,2). The above functions are also known in the small- 2: expansion [T6l[3T] - see 
Eqs. (3.9)-(3.10), (4.8) and Fig. 1 of Ref. [5]. 

The functions 0^^ (5) originate from the bremsstrahlung b -^ S'jg and b -^ S'jqq channels 
(but excluding the diagrams in Fig. [1]). Their numerical importance for b — )■ X^'y is mild due 
to the relatively high photon energy cutoff -Eq ~ ^- At large z, all the 0^^ (5) with i E {1,2} 

and j G {1,2,7,8} behave as 0{-). Consequently, they can be ignored in the next section 

(2) 
where only the leading terms of the large-2 expansion of K^j are considered. 



4 Charm decoupling: evaluation of the WCs 



-in) 



The matrices Xj," in Eq. (13. ip can be obtained by charm decoupling in the limit rric ^ ^. In 



the first step, we perform three-loop matching of the 5-fiavour theory (12. ip onto the 4-fiavour 
one given by 



c 



eff 



£ 



qcdxQedI'W? 



d,s,b) + 'Gs^v:,v,,Y.c[Q[. 



(4.1) 



Here, Q'^ differ from Qi in Eq. (12. 2p only by the absence of charm-quark currents in Q3,4,5,6- 
Additional terms containing non-physical (evanescent and/or EOM-vanishing) operators on the 
r.h.s. of Eqs. (12. ip and (14. ip are implicitly assumed. A complete list of such terms can be found 
in Sec. 3 of Ref. [H]. Here, we just quote three examples of gauge-invariant EOM-vanishing 
operators 



7sl^^^T%l)D''GI, + gQ'!^, -SLt)^^bi 
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g<'^(4.2) 



We proceed by analogy to our NNLO electroweak-scale matching for b ^ S'j and b ^ sg |13] . 
Requiring equality of appropriately renormalized off-shell Green functions in both theories leads 
to relations that allow to express Cl in terms of Cj. Expansion in external momenta is performed 
prior to loop integration. No IR regulators are introduced. All the particles except for the 
decoupled charm are treated as massless, although linear terms in rub from bottom propagators 



and operator vertices are retained. Diagrams like the one in Fig. |2] still enter our calculation, 
but going off-shell and expanding in external momenta makes their evaluation much easier. 
The necessary three-loop integrals are found with the help of MAT AD [26] . 

As well known, spurious IR divergences that appear in such a procedure (regulated dimen- 
sionally) cancel out in the final expressions for C'^ in terms of Cj. All the loop diagrams that 
contain no charm quark are scaleless, so they vanish in dimensional regularization. Thus, 1/e" 
poles on the 4-fiavour theory side originate from the UV-renormalization constants only, while 
the 5-fiavour theory poles come from loop diagrams, too. 

We identify the renormalization scale at which the matching is performed with the previously 
introduced scale fi^. However, the charm mass rric is assumed to be MS-renormalized at another 
scale called fie- The coefficients Cf^(/if,) and Cj''^^(/i{,) are expanded in terms of fig and 
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respectively, as follows: 
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cf = cr''+s^scr''+c^ci 

while Cj'*'^ are related to C'^ precisely as in Eq. (12. 3p . with the same numbers i/i and Zi. 



(4.3) 



(2) 



Once the r.h.s. of Eq. (12. 4p is perturbatively expanded up to 0{a'^), the sought K^j are 

multiplied by C| 2*7 s only. All the other unprimed WCs will be set to zero everywhere in the 
following. After such a simplification, our results for C'^'^^ take the form 
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where L/j = In-^, L^. = In-^, and ni = 3 denotes the number of fiavours that are kept 
massless throughout the calculation. Retaining Ui as a symbol is convenient for cross-checking 
the BLM-part subtraction later on. 

The above results have been obtained by matching b — )■ 57, h -^ sg and h — )■ sqq off-shell 
Green functions in both theories. As a by-product, we have also obtained WCs of EOM- 
vanishing operators like the ones in Eq. (14. 2p . However, since on-shell matrix elements of such 



6 



operators vanish [32], there is no need to consider them further. This is precisely the point 
where the Lagrangian-level decouphng is advantageous with respect to the purely diagrammatic 
approach. In the latter case, complicated on-shell integrals may occur in contributions that 
are due to EOM-vanishing operators alone, but this fact is not visible before reduction to truly 
independent master integrals. An additional advantage in our particular case is that we can 
use (in the next section) the known two-loop on-shell 6 — )■ 57 matrix element of (^4 that has 
been evaluated without reduction to master integrals [23] . 

In the remainder of this section, let us recall several important points concerning renormal- 
ization in off-shell matching calculations. First, the external fields must be renormalized in an 
identical manner on both sides of each matching equation. One possibility is to renormalize 
all the fields on shell. More conveniently, one can shift from the on-shell to the MS scheme 
on the 4-flavour theory side, and perform an identical shift on the 5-flavour side. Second, one 
adjusts the gauge coupling renormalization on the 5-flavour theory side in such a way that the 
renormalized coupling equals to a'^ exactly in e = {4 — D)/2. At one loop, the necessary 
renormalization of a'^ in the 5-fiavour theory is given by a'^^^^^ = Z'^a'^ with 



25/1 , + ,„4.^ ^r(ij^^4^Y 



Z,^l + sq--(^--7 + lntoj + -^^[-^j \ + 0(S:% (4.7) 

in full analogy to Sec. 4 of Ref. [13] where more explanations can be foundo Explicit expressions 
for shifts in the quark mass and wave function renormalization can be found there, too. 

As far as the WC renormalization is concerned ((^^^^bare _ j^jCj^Zj-), we begin with the 
MS scheme in both theories, and never redefine the Zj^ . However, we re-express them on 
the 5-fiavour side in terms of a'^ that is renormalized according to Eq. (14. 7p . This leads to 
appearance of UV-finite terms in Zij because the relation between a'^ and fig contains 0{e) 
terms. Application of D-dimensional rather than 4-dimensional relations between the gauge 
couplings has been essential for successful tests of our expressions (I4.4l) -( l4l6l) against results 
derived with the help of asymptotic expansions in the off-shell case. These tests involved direct 
three-loop 6 — ?• 57 matching between the full SM and the 4-flavour effective theory (14.11) for 
rric -C Mw, 

5 On— shell amplitudes 

We can now proceed to evaluating on-shell b -^ Xg^ amplitudes in the 4-fiavour theory using 
Cj'^*^^ as they stand in Eqs. (I4.4p -( l4^ . With all the gauge couplings factorized out and the 
overall factor of ^HfV^lVtt omitted, the relevant expressions read 
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^ Terms containing (—7 + ln47r) are often skipped in Ref. [T3]. For more details about decoupling relations 
see Ref. [33] and Eq. (12) of Ref. ^. 
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(5.2) 
(5.3) 
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where (Qj)^"'-' denotes the n-loop renormahzed matrix element of Q'- between the considered 

external states, and C^'*-'^-* are used only when C^'*-'^-* = C'^^''^'^^. 

For the penguin operators (j=3,4,5,6) we have (Qf)''^^ = {Qf)''"^-yj{Q7Y''~^^ -ZjiQ^sY"'-^^ 
with the same numbers tjj and Zj as in Eq. (12. 3p . In fact, those numbers are determined by the 



requirement {s'j\Qf\b)l.Q = = {sg\Q^j^\b)i,fl. Expressing amplitudes in terms of {Qf)es is 
convenient in other cases, too. For instance, the bremsstrahlung matrix elements {s''yg\Q^j^\b)l.Q 
are given by subsets of diagrams contributing to {s'jg\Q'j''\b)^^\ namely those where both the 
photon and the gluon are attached to the quark loop [35]. At two loops, {s^ylQ^lby^g contain 
no IR divergences, contrary to (s7|Q^''|6)^^^ 
In Eq. (15.11) . we need 
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The result in Eq. (15. 5 p has been extracted from Eqs. (3.1) and (6.21) of Ref. [23] after reintro- 
ducing explicit n^-dependence there. More precisely, setting z = in the quoted equations of 
Ref. [23] gives the same number for r^ + 747 ln(mb/ jJb) there as setting rii = A in the square 
bracket of Eq. (15. 5 p here. However, here we need n^ = 3. 



The next steps to perform are as follows: 

• Calculate moduli squared of the amplitudes in Eqs. fl5.ip - fl5.4p . sum over polarizations 
and integrate over the phase space. 

• Re-expand everything in terms of 5s using a'^ = Og (l — |5si^D + C^(5s)), and take 
into account normalization to the semileptonic rate in Eq. (12. 4p . 

• Pick up only those 0{a'^) terms that contain at least a single CI2, and read out the 
corresponding K^j . 

• Subtract the BLM contributions (leading terms in the large-2; expansion only) using 
Eqs. ( 13.2p -( l3l3l) with /3o = 11 — |(n^ + 2), and check that all the n^-terms cancel out. 

A brief look at Eqs. (I4.4p -( H3]) and (I5.ip - (l5.4p ensures that the matrix element in Eq. (15. 5p 
is actually the only two-loop on-shell one that we need. Let us stress again that this is the case 
only thanks to identifying the EOM-vanishing operators at the Lagrangian level. 

Another straightforward observation is that K^j for i,j G {1,2} receive contributions 
only from the C'j^^>^ {Q'^)'-^^ term in Eq. I^Ji) and the Ci^^^(y^'■^^ term in Eq. ([52D. As 
the latter term is not known, we shall neglect it in what followsjj With this approximation, 
everything we need is given by the quantity that multiplies (C2 — |C| j in 0!^'^^'^°^ (see 

Eq. (14. 5p ). It contains no n^-piece, and is equal to i^27 + ^(V-^)- Consequently 
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243 81 ' 
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which is identical to Eq. (5.4) of Ref. |5]. 

It remains to determine -f^mure)' ^•^•' -f^ij '^™ for i G {1,2} and j G {7,8}. Once 
(s7gg|Q4|6)*^°) has been neglected, the only relevant processes are b -^ S'j and b —^ S'jg, and 
we need at least one beyond-LO coefficient C^'^^^'^^, which makes the calculation very similar 
to the b — >■ X^'j one at the NLO. In fact, all the necessary matrix elements and phase-space 
integrals come in the same combinations as in the known results for 



-"77 


- -f^77 


= -^ + F-f^. + 40SW 


Kf^' - 


- -"78 


= f-^^^ + f^. + 2 0g(5). 


kS' - 


= 4i 


= 4</.g(5), 


-"47 


- K^S 


-2Reb{z) + ^^-lL, = Rerl 


-"48 


- kS 


= 2 0S(5). 



-|^ + >,)l, + 2 4V(5), 



(5.9) 



^ It originates precisely from the diagrams in Fig. [T] that were discussed at the end of Sec. [2] Now it is clear 
that including their n^-parts in the BLM approximation would not be mandatory because the qq pair emitted 
from Q'^ has no gluonic counterpart in any other diagram. More precisely, this counterpart occurs only in the 
first EOM-vanishing operator in Eq. ()4.2p that gives no contribution on shell. 
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For completeness, all the relevant functions (pij {5) are collected in the Appendix. Normalization 

'(1) 



to the semileptonic rate is already taken into account in K)^i' ^ and there is no other point where 



it could matter in the evaluation of -f^mvTs) • 



(2) 



Once the quantities from Eq. (15. 9p are used, equations that determine the sought Ki^ take 
a simple form 

cf^Kf^+cf^Kf^ + OiA = kSc',^'^^^ + kSc',^'^''^ + Kg^C',^'\ (5.10) 



where Cj^^^"^^ stands for C'j^'^^'^^ (14. 6 p with Cy^^^ -terms set to zero. 

In the last step, as already mentioned above, we need to subtract the BLM parts (I3.2p - (l3.3p 
from the calculated kIj to obtain K^j''^'^^. At this level, it is convenient to express K'^j' first 
in terms of K^i , and next in terms of 



^{l)rcm _ ^(1) a (26 4 r ~ 

= - W - ^ - i^^ + |Rea(l) + |Re6(l) + |f L, + ^L^ + 2 0iV(5) + O 



Our final results for K.^^vjs) ^^^^ ^^^ form 



^(2)rem _ irW irW j_ ( 127 35 r A ZV^Cl) i 2 /i r \ rA^^cm 



(5.11) 



.27 - -fV27-fV77 t y^ - 27^0 J ^'^rs "T 3 V J- - -^i;;-fV47 

4736 r 2 , 1150 r 1617980 , 20060 A , 1664 r 1 /'O 1^ 1 ^ ('c; 1 oN 

~ "729"^£' + 729"^^ ~ 19683 + ^43~^3 + "sT^^ + U (^-J , (^i).12J 

K^r^^ = K^,}^ Kil^ + {^^-§Lo)kS + 1(1 -Ln)KS + 0{^, (5.13) 

i^ff^- = -IK^^- + (^ - fi^z.) i^S - W + f Cs + il^?, - fi^D + O (i) , (5.14) 
^(2).en. ^ _l;^(2)-- + (A _ ^Ln) K^^ + O {^ , (5.15) 

which is identical to Eqs. (5.5)-(5.8) of Ref. |S]. 

6 Conclusions 

We have presented details of our large-mc calculation [5] of those NNLO corrections to 
B{B — > Xs'~f) that still require interpolation in rric. Applying Lagrangian-level decoupling 
rather than the purely diagrammatic asymptotic expansions has led to appreciable simplifica- 
tions of the analysis. 

Our results are going to be useful again in the near future when the calculation of Kij '^ 
and i^27 °™ ^t '^c = is completed [36] providing data for an upgraded interpolation in rric. 
With those inputs, as well as new results for i^7g [37] and the remaining BLM terms [20] . 
an update of the phenomenological analysis [HIS] will be mandatory. An ultimate goal is to 
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make the perturbative uncertainties in B{B — t- Xs7) 
perturbative [2] and experimental [6] ones. 



negligible with respect to the non- 



Acknowledgments 

This work has been supported by the DFG through SFB/TR9 "Computational Particle Physics" 
and the "Mercator" guest professorship programme. M.M. acknowledges partial support from 
the EU-RTN Programme "FLAVIAnet" (MRTN-CT-2006-035482), and from the Polish Min- 
istry of Science and Higher Education as a research project N N202 006334 (in years 2008-11). 
We thank the Galileo Galilei Institute for Theoretical Physics for hospitality and the INFN for 
partial support during completion of this work. 

Appendix 



Here, we quote explicit expressions for the relevant functions 



& 



^78 



'88 



(5) 



-|ln^^ 



|ln5-f + f5 + i5^ 



2^3 
9" 



\6{6-A)\n6, 



Li2(l-<5) 



1^2 
6^ 



5\n5 + \5-\5^ 



12" 



5^ + 25 + A\n{l - 5) 



2 

In^ 

mt 



+ 



4Li2(l - 5) 
-30S(5) = 



2^2 
3^ 



5(2 + (5)ln5 + 81n(l-(5) 



2rf 



35' + 7(5} . 



54 



TT 



3^3 



TT 



^^(1 



81 " 



6 + \6^ 



^^2 ^ A( 
108 ^ 54^ 



lim 



¥mi, 



-'27 



(/>) 



{^■' 



2(5 + 3) arctan 




\(S^ 



4(5 + 3 



'1-5 



arctan 



'1-5 



(A.l) 
(A.2) 

(A.3) 



(A.4) 



The functions (f)\j and 



. 3 + (5 V3 + 5 

48 are exactly the same in the 5-fiavour and 4-flavour theoriesjj They 
are generated by the s- and 6-quark loops with no Dirac traces only. Contributions with 
traces cancel out in the same way as in the QED electron-loop contributions to three-photon 
interactions (Furry theorem). 



Let us note that ^gg in Eq. (IA.3P 



i(2)/3o 
'88 5 



as well as the neglected diagrams in Fig. [T] contain 



coUinear logarithms where m^ ^ for q = u,d,s need to be retained at the perturbative 
level. The actual collinear regulators in reality are of order of the light meson masses {rriTr, 
rriK)- Non-perturbative collinear effects in 5 — >■ X^'^ have been discussed in Refs. [2l|39]. Their 
numerical effect on the branching ratio for i?o ~ -^^ is generically small (~ 1%) thanks to the 
phase-space suppression, electric charge factors and/or small values of the relevant WCs. 

^ Eq. (3.12) of Ref. [5] contains a misprint in the coefficient at Ivaim^^mt which we correct in the first fine 
of Eq. (|A.4|) here. 
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